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ABSTRACT 

Dissolved  barium  in  excess  of  the  U.S.  Environmental  Protection 
Agency  drinking  water  standard  occurs  in  groundwater  of  a  region  of  the 
Cambrian-Ordovician  aquifer  system  in  northeastern  Illinois  that  is  anoma- 
lously depleted  in  dissolved  sulfate.  The  concentration  of  dissolved 
sulfate  is  an  important  control  on  the  dissolved  barium  concentration 
because  of  the  low  solubility  of  the  mineral  barite  (barium  sulfate).  The 
depletion  of  dissolved  sulfate  is  believed  to  be  due  to  bacterial ly  mediated 
sulfate  reduction  occurring  in  zones  of  anaerobic  conditions.  Characteris- 
tics of  the  groundwater  that  support  this  interpretation  are  the  presence  of 
dissolved  hydrogen  sulfide,  increased  alkalinity,  and  31*S  enrichment  in  the 
dissolved  sulfate. 

Isotopic  analyses  determined  that  oxidation  of  sulfide  minerals  is 
the  primary  source  for  dissolved  sulfate  in  groundwater  in  the  unconfined 
region  of  the  flow  system  and  that  dissolution  of  evaporite  minerals  is  the 
primary  source  for  dissolved  sulfate  present  in  mineralized  groundwater  in 
the  confined  region  of  the  flow  system.  The  source  of  the  evaporite 
minerals  is  problematic.  Evaporites  are  not  known  to  be  present  in  the  flow 
system  and  isotopic  analyses  provided  evidence  that  the  dissolved  sulfate 
was  derived  from  evaporites  of  a  different  age  than  could  occur  in  the 
Cambrian  and  Ordovician  strata  of  the  aquifer  system. 

The  depletion  of  180  in  groundwater  in  confined  regions  of  the  flow 
system  is  interpreted  as  evidence  of  the  importance  of  glacial  processes  to 
recharge  of  the  groundwater.  Dissolved  sulfate  from  marine  evaporites  of  a 
younger  geologic  age  may  have  been  recharged  to  the  flow  system  as  a  result 
of  glacial  processes.  Further  research  is  necessary  to  establish  the  rela- 
tion between  glacial  events  and  the  chemical  composition  of  groundwater  in 
confined  regions  of  the  Cambrian-Ordovician  aquifer  system. 

KEY  WORDS 

aquifer  properties/barite/barium/chemical  properties/evaporite  minerals/ 
glacial  events/groundwater/microbial  reactions/sedimentary  rocks/stable 
isotopes/sulfate 
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INTRODUCTION 

Chemical  analyses  of  groundwater  samples  performed  by  the  Illinois 
Environmental  Protection  Agency  determined  that  the  dissolved  barium  concen- 
tration in  a  large  number  of  wells  in  Illinois  exceeds  the  drinking  water 
standard  of  1.0  mg/L  (U.S.  EPA,  1973).  The  majority  of  wells  exceeding  the 
standard  are  located  in  northeastern  Illinois.  Research  at  the  Illinois 
State  Geological  Survey  (Gilkeson  et  al . ,  1978)  established  (1)  that  the 
affected  wells  in  northern  Illinois  are  all  finished  in  bedrock  of  the  Cam- 
brian and  Ordovician  Systems,  and  (2)  that  the  concentration  of  sulfate 
dissolved  in  groundwater  is  an  important  control  on  the  concentration  of 
dissolved  barium. 

For  a  significant  number  of  groundwater  samples,  the  ion-activity 
product  of  barium  (Ba2+)  and  sulfate  (SC\2~)  was  in  equilibrium  with  the 
mineral  phase  barite  (BaSC\)  according  to  the  following  reaction: 

BaSC\  *fc  Ba2+  +  SC\2-  (1) 

This  reaction  requires  an  inverse  relation  between  the  activities  of  Ba2+ 
and  SC\2"  ions.  High  barium  concentrations  may  occur  only  in  groundwater 
that  has  a  low  concentration  of  dissolved  sulfate.  In  a  large  region  of 
northeastern  Illinois,  high  barium  concentrations  occur  in  groundwater  that 
is  anomalously  depleted  in  dissolved  sulfate. 
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The  depleted  groundwater  also  commonly  contains  dissolved  hydrogen 
sulfide.  The  reaction  of  sulfate  ion  (SO,2-)  with  organic  matter  (CH20-R), 
where  R  is  an  undefined  radical,  to  form  bicarbonate  ion  (HC03")  and  hydro- 
gen sulfide  (H2S)  can  be  written  as  follows: 

2(CH20-R)  +  SO,2"  -+■   2(HC03-)  +  H2S  +  R  (2) 

Although  thermodynamic  calculations  predict  that  the  inorganic 
reduction  of  sulfate  to  sulfide  is  possible  at  the  temperatures  and  pres- 
sures that  exist  in  the  aquifer  system,  it  is  well  known  from  laboratory 
experiments  that  the  inorganic  reaction  will  not  occur  in  the  aquifer 
system.  The  inorganic  reaction  has  been  accomplished  in  the  laboratory 
only  at  temperatures  greater  than  250°C  (Goldhaber  and  Kaplan,  1974). 
Given  the  conditions  in  the  aquifer  system,  however,  the  reduction  of  sul- 
fate has  been  known  to  occur  only  through  bacterial  mediation  of  the  reac- 
tion. Sul fate-reducing  bacteria  require  an  anaerobic  environment  and  a 
source  of  organic  carbon  (Stumm  and  Morgan,  1970).  The  chemical  composi- 
tion of  groundwater  from  rocks  of  the  Cambrian  and  Ordovician  Systems  of 
northern  Illinois  was  studied  to  determine  if  the  depletion  of  dissolved 
sulfate  could  be  the  result  of  bacterial  reduction.  Stable  isotopes  of 
sulfur  were  analyzed  because  the  bacterial  reduction  of  sulfate  results  in 
a  large  fractionation  of  sulfur  isotopes  between  sulfate  and  sulfide  (Clay- 
pool  et  al. ,  1980). 

The  oxidation-reduction  potential  (Eh),  the  titration  alkalinity, 
and  the  concentration  of  dissolved  organic  carbon  in  the  groundwater  were 
also  studied.  An  attempt  was  made  to  culture  anaerobic  bacteria  in  ground- 
water from  selected  well  sites,  and  stable  isotopes  of  hydrogen  and  oxygen 
in  the  water  molecule  were  analyzed  to  provide  information  on  rock-water 
interactions  and  the  history  of  groundwater  in  the  flow  system. 

HYDROGEOLOGIC  SETTING 

The  location  of  the  study  area  in  northeastern  Illinois  is  shown  in 
figure  1.  A  generalized  column  of  rock-stratigraphic  units  and  aquifers  in 
northeastern  Illinois  is  shown  in  figure  2.  The  terminology  defining 
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FIGURE  1.  Map  of  the  study  area  showing  distribution  of  the  Maquoketa  Shale  Group,  approxi- 
mate direction  of  groundwater  flow  in  the  Cambrian-Ordovician  aquifer  system  before  large  ground- 
water withdrawals,  zones  of  strong  reducing  potential,  and  well  sites  discussed  in  the  text. 
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Generalized  column  of  rock  stratigraphic  units  and  aquifers  in  northeastern 
(Prepared  by  M.   L.   Sargent,    Illinois  State  Geological   Survey,   1981.) 


aquifer  systems  is  from  Hughes,  Kraatz,  and  Landon  (1969).  The  sequence 
from  the  Galena-Platteville  Dolomite  through  the  Ironton-Galesville  Sand- 
stone functions  hydraulically  as  one  aquifer  (Walton,  1965)  and  is  referred 
to  as  the  Cambrian-Ordovician  aquifer  system.  Depending  upon  local  condi- 
tions and  groundwater  needs,  wells  penetrate  the  aquifer  system  to  differ- 
ent depths.  The  main  water-yielding  units  are  the  Glenwood-St.  Peter  Sand- 
stone and  the  Ironton-Galesville  Sandstone. 

Shales  in  the  Maquoketa  Shale  Group  and  the  Eau  Claire  Formation  are 
important  confining  units.  The  shales  in  the  Eau  Claire  Formation  separate 
the  Cambrian-Ordovician  aquifer  system  from  deeper  water-bearing  units. 
The  Elmhurst  Sandstone  Member  at  the  base  of  the  Eau  Claire  Formation  and 
the  Mt.  Simon  Sandstone  are  connected  hydraulically  and  form  the  Mt.  Simon 
aquifer  system.  No  wells  receive  groundwater  exclusively  from  the  deeper 
aquifer  system;  wells  usually  penetrate  only  a  few  hundred  feet  into  this 
aquifer  because  the  water  quality  decreases  with  depth. 

In  north-central  Illinois,  where  the  Maquoketa  Shale  Group  is  eroded, 
recharge  of  the  Cambrian-Ordovician  aquifer  system  is  produced  from  verti- 
cal leakage  of  meteoric  water  through  glacial  drift.  The  regional  distri- 
bution of  the  Maquoketa  Shale  Group  is  shown  in  figure  1.  In  northeastern 
Illinois,  where  the  Maquoketa  strata  overlies  and  confines  the  aquifer, 
much  recharge  is  produced  by  vertical  leakage  through  the  shale  from  rocks 
of  Silurian  age  or  the  glacial  drift  (Walton,  1965). 

The  exact  configuration  of  the  potentiometric  surface  of  the  Cam- 
brian-Ordovician aquifer  system  before  extensive  groundwater  withdrawal  is 
unknown.  Interpretation  of  early  data  indicates  that  under  natural  condi- 
tions the  surface  was  relatively  featureless  and  sloped  gently  to  the 
southeast  (Suter  et  al . ,  1959).  Significant  areas  of  recharge  were  in 
Boone  and  De  Kalb  Counties.  Because  of  their  topographic  relief,  major 
drainageways  such  as  the  Fox,  Des  Plaines,  and  Illinois  Rivers  must  have 
had  a  significant  effect  on  hydraulic  gradients  of  the  natural  potentiomet- 
ric surface. 

As  a  result  of  large  withdrawals  of  groundwater,  broad,  deep  cones 
of  depression  have  formed  in  Cook,  Du  Page,  and  Will  Counties.  The  steep- 
ened hydraulic  gradients  have  increased  the  transmittal  of  groundwater 


toward  the  centers  of  heavy  withdrawal  (Sasman  et  al . ,  1977).  Groundwater 
flow  is  from  all  directions  toward  these  pumping  centers. 

In  the  Chicago  region,  the  potentiometric  surface  of  the  aquifer  has 
changed  more  than  210  meters  (700  ft)  in  some  localities  (Sasman  et  al . , 
1977).  Significant  dewatering  of  the  Galena-Platteville  Dolomite  has 
occurred  in  Will,  Cook,  Du  Page,  and  eastern  Kane  Counties. 

HYDROGEOCHEMISTRY 

The  location  of  the  study  area  is  shown  in  figure  1.  Data  on  the 
chemistry  of  groundwater  at  well  sites  shown  on  figure  1  is  presented  in 
table  1A.  Data  on  well  construction  is  presented  in  table  IB.  In  the 
region  where  the  Maquoketa  Shale  Group  is  eroded  and  significant  recharge 
occurs  from  infiltration  of  meteoric  waters,  the  major  ions  in  groundwater 
are  calcium,  magnesium,  and  bicarbonate.  Eastward  and  southeastward  from 
this  region  the  groundwater  becomes  more  mineralized  along  the  probable 
direction  of  predevelopment  flow. 

An  estimated  predevelopment  flow  line  is  shown  by  the  shaded  arrow 
in  figure  1.  The  flow  line  is  based  on  an  interpretation  of  the  predevel- 
opment potentiometric  surface  (Suter  et  al . ,  1959)  and  does  not  consider 
the  effects  that  anisotropy  or  heterogeneity  in  the  rock  strata  may  have  on 
the  direction  of  groundwater  flow.  The  flow  line  is  intended  only  as  an 
approximation  of  groundwater  flow  in  the  aquifer  system  to  provide  a  frame- 
work for  the  study  of  regional  variation  in  the  composition  of  groundwater. 
The  locations  of  well  sites  along  the  flow  line  are  shown  on  traverse  A-A1 
in  figure  1.  Data  on  the  chemical  composition  of  the  groundwater  along 
the  traverse  A-A'  are  presented  in  table  2  and  figure  3.  The  distribution 
of  dissolved  sulfate  in  groundwater  of  the  Cambrian-Ordovician  aquifer  sys- 
tem for  a  region  of  northeastern  Illinois  is  shown  in  figure  4.  Over  a 
large  area  of  northeastern  Illinois,  the  dissolved  sulfate  concentration  in 
the  aquifer  system  is  less  than  10  mg/L.  East  and  southeastward  of  this 
region,  the  dissolved  sulfate  concentration  increases  rapidly  and  is  pre- 
sent at  concentrations  of  greater  than  800  mg/L  in  southeastern  Cook  County. 
Sodium  and  chloride  are  other  chemical  constituents  that  are  present  in 
significant  concentrations  as  the  groundwater  becomes  more  mineralized. 
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FIGURE  4.  Distribution  of  dissolved  sulfate  in  groundwater  of  the  Cambrian-Ordovician  aquifer 
system  for  a  region  of  northeastern  Illinois.  (From  Gilkeson  and  Cartwright,  in  press.) 
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The  distribution  of  dissolved  barium  in  groundwater  of  the  Cambrian- 
Ordovician  aquifer  system  for  a  region  of  northeastern  Illinois  is  shown  in 
figure  5,  and  the  variation  in  dissolved  barium  that  occurs  along  flow  line 
A-A'  is  shown  in  figure  3.  Groundwater  from  the  Cambrian-Ordovician  aqui- 
fer system  exceeds  the  drinking  water  standard  of  1.0  mg/L  for  barium  over 
a  large  part  of  McHenry  and  Kane  Counties,  the  western  part  of  Lake  County, 
the  northwestern  part  of  Du  Page  County,  and  the  northeastern  part  of  De 
Kalb  County.  At  some  well  sites  the  dissolved  barium  concentration  exceeds 
15  mg/L.  Groundwater  with  barium  concentrations  greater  than  1.0  mg/L  gen- 
erally have  dissolved  sulfate  concentrations  less  than  10  mg/L. 

A  common  observation  in  confined  groundwater  flow  systems  is  a  gen- 
eral decrease  in  the  oxidation-reduction  (redox)  potential  of  groundwater 
in  the  direction  of  flow.  The  decline  in  redox  potential  is  believed  to  be 
caused  by  a  sequence  of  oxidation-reduction  reactions  that  reduces  oxidized 
species  initially  present  in  the  groundwater.  These  reactions  are  dis- 
cussed by  Stumm  and  Morgan  (1970);  Champ,  Gulens,  and  Jackson  (1979);  and 
Cherry  et  al .  (1979). 

To  characterize  the  variation  in  redox  potential  of  groundwater  in 
the  Cambrian-Ordovician  aquifer  system,  potentials  were  measured  using  a 
platinum  electrode  coupled  to  a  calomel  reference  electrode  mounted  in  an 
airtight  chamber.  The  measurement  technique  is  discussed  by  Back  and 
Barnes  (1965).  In  the  present  study,  potentials  were  monitored  using  two 
platinum  electrodes  to  guard  against  electrode  "poisoning."  Redox  poten- 
tials of  groundwater  presented  in  table  1A  are  expressed  as  Eh  in  milli- 
volts in  terms  of  the  normal  hydrogen  electrode.  The  platinum  electrode 
values  are  used  only  in  a  qualitative  sense  to  evaluate  trends  that  are 
present  within  the  flow  system;  the  values  are  not  to  be  interpreted  within 
the  framework  of  equilibrium  thermodynamics. 

In  the  western  part  of  the  study  area,  Eh  values  are  relatively  high 
where  oxidized  meteoric  water  recharges  the  Cambrian-Ordovician  aquifer 
system  through  a  thin  cover  of  glacial  drift.  Examples  are  Lee-1,  Ogle-2, 
Ogle-3,  and  Boone-1.  As  groundwater  moves  southeastward  and  becomes  con- 
fined by  fine-grained  glacial  drift  and  the  Maquoketa  Shale  Group,  there  is 
a  decrease  in  the  redox  potential.  The  change  in  redox  potential  that 
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FIGURE  5.  Distribution  of  dissolved  barium  in  groundwater  of  the  Cambrian-Ordovician  aquifer 
system  for  a  region  of  northeastern  Illinois.   (From  Gilkeson  and  Cartwright,  in  press.) 
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occurs  along  flow  path  A-A'  (fig.  1)  is  shown  on  the  traverse  A-A'  in 
figure  3.  The  redox  decreases  with  the  direction  of  groundwater  flow  to  a 
zone  of  strong  reducing  potential  in  Kane  County.  The  cross-hatched  pattern 
in  figure  1  delimits  the  region  known  to  be  a  strong  reducing  zone.  Eh 
values  as  low  as  -100  millivolts  have  been  measured  for  groundwater  within 
this  zone.  From  central  Kane  County  the  zone  extends  northward  along  the 
Fox  River  Valley  into  McHenry  County.  The  extent  of  this  region  in  McHenry 
County  cannot  be  accurately  determined  from  the  available  information. 

Anaerobic  conditions  also  exist  in  localized  areas  where  the 
Maquoketa  is  eroded  and  the  flow  system  is  confined  by  fine-grained  glacial 
drift.  Examples  are  well  sites  De  Kalb-1  and  De  Kalb-2  (fig.  1  and 
table  1A).  The  regional  extent  of  anaerobic  conditions  at  these  well  sites 
is  not  well  known. 

One  indication  that  bacterial ly  mediated  sulfate  reduction  may  be 
occurring  in  the  anaerobic  zones  is  the  presence  of  hydrogen  sulfide  in  the 
groundwater.  The  reduction  of  sulfate  (S042~)  results  in  the  generation  of 
hydrogen  sulfide  (H2S)  and  bicarbonate  (HC03~)  as  demonstrated  in  equa- 
tion 2.  In  the  range  of  pH  present  in  the  groundwater  of  the  aquifer  sys- 
tem, a  significant  amount  of  the  hydrogen  sulfide  (H2S)  dissociates  to  form 
bisulfide  (HS").  The  major  source  of  titration  alkalinity  in  the  aquifer 
system  is  bicarbonate  (HC03~)  from  the  dissolution  of  carbonate  minerals. 
Bacterial ly  mediated  sulfate  reduction  would  be  expected  to  increase  the 
titration  alkalinity  in  groundwater  through  addition  of  HC03"  and  HS". 
Thus  the  ratio  of  titration  alkalinity  to  hardness  would  increase  and  act 
as  an  indicator  of  the  process.  The  ratio  of  titration  alkalinity  to  hard- 
ness (as  equivalents  of  calcium  carbonate)  is  shown  for  groundwater  samples 
along  traverse  A-A'  in  figure  3.  The  titration  alkalinity  is  enhanced  rela- 
tive to  hardness  in  groundwater  samples  from  the  reducing  zone  and  the  tran- 
sition zone  downgradient  of  the  reducing  zone.  As  groundwater  becomes  more 
mineralized  along  the  flow  path  there  is  a  decline  in  the  ratio  of  titration 
alkalinity  to  hardness  (well  site  C-4  to  C-6,  fig.  3).  The  decline  is  pos- 
sibly due  to  precipitation  of  carbonate  minerals  in  this  region  of  the  flow 
system.  Carbonate  mineral  equilibrium  reactions  in  the  flow  system  have 
not  been  investigated  but  are  expected  to  be  studied  in  the  future. 
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Analyses  of  total  organic  carbon  in  groundwater  samples  from  well 
sites  in  the  study  area  are  shown  in  table  1A.  Organic  carbon  is  generally 
present  in  samples  from  well  sites  located  in  the  unconfined  region  of  the 
aquifer  system  in  north-central  Illinois.  Concentrations  as  great  as 
21.0  mg/L  have  been  measured  at  well  site  Lee-1.  The  concentration  of 
total  organic  carbon  decreases  along  flow  paths,  and  the  measured  concentra- 
tions from  well  sites  located  in  the  strong  reducing  zone  were  0.7  mg/L  at 
well  site  Kane-2  and  0.8  mg/L  at  well  site  Kane-3.  Well  sites  located  in  a 
region  of  oxidized  groundwater  downgradient  from  the  strong  reducing  zone 
contained  total  organic  carbon  in  concentrations  less  than  the  limit  of 
detection  (<0.1  mg/L). 

Experimental  research  conducted  by  Zobell  and  Grant  (1942)  determined 
that  anaerobic  microbial  requirements  of  organic  carbon  are  less  than 
0.1  mg/L.  Analyses  performed  in  the  present  research  indicate  that  ground- 
water in  the  reducing  zone  has  a  sufficient  concentration  of  organic  carbon 
to  support  microbial  activity.  Attempts  to  culture  sulfate-reducing  bac- 
teria in  groundwater  from  wells  within  the  strong  reducing  zone  have  been 
unsuccessful  to  date.  However,  methanogens  were  cultured  in  groundwater 
collected  from  well  site  Kane-2  (Risatti,  1981,  Illinois  State  Geological 
Survey,  personal  communication).  Groundwater  from  many  of  the  wells  in  the 
strong  reducing  zone  contain  slight  amounts  of  methane.  Isotopic  analyses 
of  carbon  in  the  methane  determined  that  the  methane  was  generated  by  micro- 
bially  mediated  reactions  under  anaerobic  conditions  (Coleman,  1981,  Illi- 
nois State  Geological  Survey,  personal  communication). 

STABLE  ISOTOPE  GEOCHEMISTRY 

During  this  study,  isotope  ratios  were  determined  for  the  water  and 
the  sulfate  molecules:  The  ratio  of  oxygen  -18  to  oxygen  -16  [6180/u  n\]  and 


the  ratio  of  hydrogen  -2  (deuterium)  to  hydrogen  -1  [<5D,„  QJ  were  calculated 
for  the  water  molecule.  The  ratio  of  sulfur  -34  to  sulfur  -32  [<531+S/SQ  2-\] 
and  the  ratio  of  oxygen  -18  to  oxygen  -16  [6180/<~n  2-\]   were  calculated  for 
the  sulfate  molecule.  The  isotope  ratios  are  expressed  in  delta  units  (6)  as 
per  mille  (parts  per  thousand,  °L)   differences  related  to  a  standard.  In 
this  research,  the  <5D  and  6180  are  expressed  as  per  mille  differences 
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compared  to  the  VSMOW  standard  of  ocean  water  (Gonfiantini ,  1978),  and  the 
631*S  analyses  referenced  to  Canyon  Diablo  meteoritic  troilite  (Rees,  1978). 
The  equation  for  calculating  the  isotope  ratio  is 

"•  "  [(Rsamp  "  "W  '  Rstd'  *  100°  (3) 

where  RsamD  and  R  ..  are  the  isotope  ratios  D/H,  180/i6q  ,  *kS/32$    of  the 
sample  and  the  standard  respectively. 

Information  on  the  sampling  procedure  and  sample  preservation  is 
presented  in  a  paper  by  Perry,  Grundl ,  and  Gilkeson  (in  press).  A  detailed 
discussion  of  analytical  procedures  is  presented  in  a  Master's  thesis  by 
Grundl  (1980). 

IS0T0PIC  VARIATION  IN  THE  WATER  MOLECULE 

Evaporation,  condensation,  and  precipitation  processes  have  a  very 
important  effect  on  the  isotope  ratio  of  meteoric  water  (Dansgaard,  1964). 
Also,  because  both  condensation  and  isotope  fractionation  are  strongly 
affected  by  temperature,  the  isotopic  composition  of  precipitation  is 
temperature  dependent. 

Water  vapor  from  evaporation  of  the  ocean  is  depleted  in  180  by 
approximately  12%  to  15%  and  in  D  by  80%  to  120%  relative  to  ocean 
water.  When  the  water  vapor  condenses,  the  rain  or  snow  that  forms  is  en- 
riched in  the  heavy  isotopes  relative  to  the  remaining  water  vapor.  As  the 
water  vapor  is  circulated  in  northward-flowing  continental  atmospheric  sys- 
tems, the  recurrence  of  evaporation-condensation  events  results  in  a  lati- 
tudinal trend  in  the  stable  isotope  composition  of  the  water  molecule.  The 
6D  and  6180  composition  of  precipitation  correlates  on  a  global  basis 
according  to  the  equation 

SD%  =  -8  6180  +  10  (4) 

that  is  known  as  the  meteoric  water  line  (Friedman,  1953;  Craig,  1961; 
Dansgaard,  1964).  For  further  information  on  processes  that  isotopically 
fractionate  meteoric  water,  the  reader  is  referred  to  a  paper  by  Dansgaard 
(1964). 

The  isotopic  composition  of  meteoric  water  is  preserved  after  the 
water  infiltrates  the  soil  and  enters  a  groundwater  flow  system.  At  the 
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high  temperatures  (50°C  to  100°C)  present  in  deep  groundwater  flow  systems, 
the  isotopic  composition  of  water  may  be  altered  through  chemical  reaction 
with  the  aquifer  rock. 

6D  and  618  values  for  brines  in  the  Illinois,  Michigan,  and  Alberta 
basins  were  determined  during  research  by  Clayton  et  al .  (1966).  Ground- 
water samples  in  each  basin  were  found  to  have  a  large  range  of  180/160 
ratios.  However,  for  each  basin,  the  least  saline  groundwaters  approached 
the  meteoric  water  line;  the  more  saline  groundwaters  departed  from  the 
meteoric  water  line  in  the  direction  of  180  enrichment.  The  high  tempera- 
tures in  the  brines  (temperatures  as  high  as  60°C  in  the  Illinois  Basin) 
were  believed  to  be  the  major  control  on  the  180  enrichment  through  exchange 
with  the  aquifer  strata. 

Table  3  presents  analyses  for  6D/H  ~s   and  <5180/n  q\  for  groundwater 
samples  collected  from  eight  well  sites  in  Illinois  for  this  research. 
Water  temperatures  at  the  well  head  ranged  from  12.0°C  to  22.8°C.  Total 
dissolved  minerals  as  filterable  residue  ranged  from  269  mg/L  to  2520  mg/L. 
The  range  in  temperature  and  total  dissolved  minerals  between  the  eight 
well  sites  encompasses  the  variation  that  occurs  in  the  Cambrian-Ordovician 
aquifer  system  in  Illinois  where  it  is  used  as  a  source  of  potable  water. 
The  close  relationship  of  the  isotope  data  from  the  well  sites  to  the  mete- 
oric water  line  is  shown  in  figure  6.  Because  of  the  distribution  of  the 
data  on  the  meteoric  water  line,  Perry,  Grundl ,  and  Gilkeson  (in  press)  con- 
cluded that  exchange  between  the  water  molecule  and  the  aquifer  rock  strata 
has  been  insignificant  and  that  groundwater  in  the  potable  region  of  the 
aquifer  system  has  retained  a  meteoric  isotopic  composition.  The  variation 
in  6D/H  Qv  and  <5180/H  n\  that  occurs  in  the  flow  system  is  concluded  to  be 
caused  by  the  variation  in  climatic  conditions  occurring  when  the  water  was 
recharged  to  the  flow  system.  Contemporary  meteoric  water  in  northern  Illi- 
nois has  an  average  <5180/H  0)  of  aDOut  -7%°  (Graf,  Friedman,  and  Meents, 
1965).  Groundwater  samples  collected  for  this  study  that  were  found  to  be 
greatly  depleted  in  180  compared  to  the  contemporary  value,  were  recharged 
to  the  flow  system  during  a  much  colder  climate.  The  most  depleted  samples 
are  from  the  southwestern  portion  of  the  aquifer  system  (Peoria,  Fulton, 
McDonough,  and  Hancock  Counties).  6180  values  of  -11%0  to  -12%  are  common 
in  this  region. 
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FIGURE  6.  Covariation  of  180  and  D  in  groundwater  from  selected  well  sites  in  northern  Illinois. 
The  close  relationship  of  the  data  to  the  meteoric  water  line  is  shown.  Data  on  individual  well 
sites  are  presented  in  table  3.  (Modified  from  Perry,  Grundl ,  and  Gilkeson,  in  press.) 
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The  variation  in  <$180/H  q)  tnat  occurs  along  a  flow  line  (A-A1)  of 
the  pre-development  flow  system  in  northeastern  Illinois  is  shown  in  fig- 
ure 3.  The  values  range  from  -7.24%0  where  the  aquifer  system  is  uncon- 
fined  in  north-central  Illinois  (well  site  Boone-1)  to  -9.67%0  in  south- 
eastern Cook  County  (well  site  Cook-6).  The  isotopic  data  indicate  that 
groundwater  in  the  unconfined  region  of  the  aquifer  system  in  north-central 
Illinois  has  an  oxygen  isotopic  composition  nearly  identical  to  contemporary 
meteoric  water.  In  the  confined  region  of  the  aquifer  system,  the  depletion 
of  180  is  compatible  with  an  increased  residence  time  of  groundwater. 
Groundwater  that  is  significantly  depleted  in  180  requires  a  source  of 
recharge  related  to  glaciation. 

Problems  with  determining  the  source  of  the  isotopically  depleted 
groundwater  in  the  flow  system  and  establishing  a  temporal  framework  for 
the  groundwater  on  the  basis  of  its  180  isotopic  composition  are  summarized 
in  a  discussion  by  Perry,  Grundl ,  and  Gilkeson  (in  press).  During  the 
Pleistocene  Epoch  the  sources  of  recharge  to  the  aquifer  system  were  mete- 
oric water  and  glacial  melt  water.  The  relative  contribution  of  each  source 
is  not  known  and  probably  varied  over  time  as  a  function  of  the  glacial  ice 
regimen.  The  isotopic  composition  of  meteoric  water  in  atmospheric  circula- 
tion systems  over  the  Midwest  during  the  Pleistocene  is  not  well  understood. 
The  actual  isotopic  composition  of  the  glacial  ice  is  also  unknown. 

It  is  also  difficult  to  determine  a  temporal  framework  for  the  ground- 
water based  on  its  oxygen  isotopic  composition.  A  common  conclusion  from  a 
variety  of  studies  (ice  cores,  ocean  bottom  sediments,  and  faunal  assem- 
blages) is  that  the  Wisconsinan  Stage  of  the  Pleistocene  Epoch  ended  \/ery 
rapidly  within  a  3,000-year  interval  at  about  13,000  years  before  the  pre- 
sent. The  present  Midwest  climate  regime  may  have  become  established  10,000 
years  ago.  In  north-central  Illinois,  where  the  aquifer  system  is  uncon- 
fined, the  residence  time  of  groundwater  with  contemporary  isotopic  composi- 
tion may  range  from  contemporary  to  several  thousand  years,  possibly  as 
great  as  10,000  years.  Groundwater  in  the  confined  region  of  the  aquifer 
system  that  is  significantly  depleted  in  180  may  have  a  residence  time 
greater  than  13,000  years. 
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ISOTOPIC  VARIATION  OF  THE  DISSOLVED  SULFATE  MOLECULE 

In  natural  systems,  sulfur  is  cycled  between  sulfate  and  sulfide 
reservoirs  (Holser  and  Kaplan,  1966).  Given  the  conditions  of  temperature 
and  pressure  in  the  Cambrian-Ordovician  aquifer  system  in  northeastern  Illi- 
nois, the  reduction  of  sulfate  to  sulfide  requires  anaerobic  biologic  medi- 
ation. During  reduction  (equation  2),  the  light  sulfur  isotope  (32S)  is 
preferentially  reduced  to  sulfide.  The  result  is  that  the  remaining  dis- 
solved sulfate  reservoir  becomes  enriched  in  the  heavier  isotope  3I+S. 
Claypool  et  al .  (1980)  concluded  from  a  review  of  research  on  bacterial  re- 
duction of  sulfate  that  a  fractionation  between  the  sulfate  and  the  sulfide 
phase  of  approximately  40%o  is  common  in  natural  systems.  Large  fractiona- 
tion effects  of  sulfur  isotopes  are  known  only  to  occur  during  the  bacterial 
reduction  of  sulfate  to  sulfide.  Significant  fractionation  of  sulfur  iso- 
topes does  not  occur  during  the  nonbiological  reduction  of  sulfate  to  sul- 
fide (Goldhaber  and  Kaplan,  1974)  or  during  the  dissolution  or  precipita- 
tion of  sulfur-bearing  minerals  (Claypool  et  al.,  1980). 

Table  1A  presents  analyses  for  631*S/™  2-\   in  groundwater  samples 
from  the  Cambrian-Ordovician  aquifer  system  in  northeastern  Illinois.  For 
well  sites  in  the  study  area,  the  <53I+S/S0  z-\   of  oxidized  groundwater  in 
the  unconfined  region  of  the  aquifer  system  ranges  from  +2.40%  to  +10.34%. 
Analyses  of  samples  from  seven  well  sites  in  the  oxidized,  unconfined 
region  of  the  aquifer  system  in  north-central  Illinois  were  used  to  deter- 
mine a  median  631tS/SQ  2-\  of  5.07%  with  a  range  of  2.36%  to  10.34%  (data 
from  this  research  and  Perry,  Grundl ,  and  Gilkeson  [in  press]).  The  sul- 
fate in  oxidized  groundwater  in  the  unconfined  flow  system  is  isotopically 
light;  and  Perry,  Grundl,  and  Gilkeson  (in  press)  concluded  that  an  impor- 
tant source  of  the  sulfate  was  the  oxidation  of  sulfide  minerals  dissolved 
during  the  recharge  of  meteoric  water. 

An  additional  minor  source  of  sulfate  is  from  precipitation  recharged 
to  the  aquifer  system.  The  concentration  and  isotopic  composition  of  sul- 
fate in  rainfall  over  a  region  of  the  study  area  (southern  Du  Page  County) 
were  investigated  by  Holt,  Engelkeimer,  and  Venters  (1972).  The  average 
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63**S/<.n  2-\  of  dissolved  sulfate  in  rainfall  was  +5%;  the  concentration  of 
dissolved  sulfate  ranged  from  1  to  10  mg/L.  The  results  are  similar  to  the 
observation  of  Rightmire  et  al .  (1974)  that  dissolved  sulfate  concentrations 
in  rainfall  over  continental  areas  usually  range  from  1  to  10  mg/L,  with  an 
average  of  3  mg/L.  The  available  data  indicate  that  the  dissolved  sulfate 
derived  from  precipitation  has  a  <53I+S/S0  2-\  that  is  nearly  identical  to 
that  of  dissolved  sulfate  derived  from  dissolution  and  oxidation  of  sulfide 
minerals. 

The  variation  in  <531*S/S0  2-\  that  occurs  along  the  flow  line  A-A'  in 
the  predevelopment  flow  system  in  northeastern  Illinois  is  shown  in  figure  3. 
A  large  fractionation  of  <53ltS/so  2  .  is  apparent  between  oxidized  ground- 
water in  the  unconfined  region  of  the  flow  system  and  groundwater  samples 
from  well  sites  located  in  the  strong  reducing  zone  where  the  aquifer  sys- 
tem is  confined  by  Maquoketa  Shale.  The  dissolved  sulfate  in  the  reducing 
zone  is  enriched  in  3kS   relative  to  the  average  composition  of  dissolved 
sulfate  in  the  upgradient,  unconfined  region  of  the  aquifer  system  with 
634S/S0  2-)  values  of  +40.30%  at  well  site  Kane-2  and  +34.11%  at  well  site 
Kane-3.  The  significance  of  the  enrichment  with  regard  to  bacterial ly  medi- 
ated sulfate  reduction  is  discussed  later  in  this  paper. 

In  groundwater  that  is  downgradient  of  the  reduced  zone,  there  is  a 
decrease  in  the  <$31*S/S0  2-\  corresponding  to  a  marked  increase  in  both  the 
dissolved  sulfate  concentration  and  the  oxidation-reduction  potential.  How- 
ever, the  <53I+S/S0  2-\  of  dissolved  sulfate  in  this  region  of  the  aquifer 
system  is  too  enriched  in  3t*S  (+23.40%  at  Du  Page-2  and  +19.56%  at  Cook-6) 
to  have  been  derived  from  oxidation  of  sulfide  minerals.  Perry,  Grundl ,  and 
Gilkeson  (in  press)  concluded  that  a  major  portion  of  the  dissolved  sulfate 
in  groundwater  of  this  region  was  derived  from  dissolution  of  evaporite  min- 
erals, but  recognized  that  the  source  of  the  evaporites  was  problematic. 
The  sulfur  isotopic  composition  of  marine  sulfates  has  varied  systematically 
throughout  geologic  time  as  a  function  of  the  distribution  of  sulfur  in  the 
sulfate  and  sulfide  reservoirs  (Holser  and  Kaplan,  1966.  Figure  7  presents 
631*S/^q  2-\  values  from  worldwide  analyses  of  evaporite  sulfates  of  differ- 
ent geological  ages.  Average  <531*S/C;n  2-\  values  for  marine  evaporites  of 
North  America  are  approximately  +30%  in  the  Cambrian  and  Ordovician  Systems 
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FIGURE  7.  Summary  of  sulfur  isotope-age  curve  for  sulfate.  Solid  areas  or  lines 
qualitatively  indicate  the  number  of  analyses,  plotted  at  their  most  probable  age. 
Horizontal  dashed  lines  signify  the  range  of  relatively  few  analyses  (Claypool 
et  al.,  1980). 
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5      S(so42-)  =  1-26  5  180(so42-)  +0.i 
correlation  coefficient,  r2  =  0.94 


0-3  0    Well  sites  with  oxidized  ground  water.  Data  in  table  1 A 

Oxidized  ground  water  in  northwestern  Illinois. 
Data  from  Perry,  Grundl,  and  Gilkeson  (in  press) 

C-2  O    Well  sites  with  hydrogen  sulfide.  Data  in  table  1 A 
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FIGURE  8.  Covariation  of  631*S  and  6180  in  dissolved  sulfate  in  groundwater  of  the  Cambrian- 
Ordovician  aquifer  system  in  northern  Illinois.  The  linear  relationship  is  for  sulfide-free 
groundwater  samples.  The  four  sulfide-bearing  groundwater  samples  that  are  displaced  from  the 
line  are  interpreted  to  be  enriched  in  3US  due  to  bacterial  sulfate  reduction. 
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(Claypool  et  al.,  1980).  The  <S34S/SQ  2-\  of  sulfi de-free  groundwater  sam- 
ples in  the  flow  system  are  not  sufficiently  enriched  in  3l|S  to  indicate 
that  sulfur  is  obtained  from  marine  evaporites  of  Cambrian  or  Ordovician 
age.  An  additional  problem  is  that  marine  evaporites  are  not  known  to  be 
present  in  the  Cambrian  and  Ordovician  strata  in  northern  Illinois.  Another 
source  of  this  dissolved  sulfate  is  proposed  in  the  discussion  section  of 
this  paper. 

The  covariation  of  <534S/S0  2-\  and  <5180/so  2-\  in  groundwater  sam- 
ples of  the  aquifer  system  is  shown  in  figure  8.  The  relationship  between 
63US/S0  2-\   and  <$180/SQ  2-\  for  the  sulfide-free  groundwater  samples  is 
described  by  the  relation 

6"S(S0^)  ■  l-M«»0(SOtl.,  ♦  0.69  (5) 

with  a  correlation  coefficient,  r2 ,  of  0.94.  The  relationship  includes  sul- 
fate derived  from  oxidation  of  sulfide  minerals  and  sulfate  derived  from 
dissolution  of  sulfate  minerals.  Although  the  mechanism  responsible  for 
the  observed  relationship  between  the  sulfide-free  groundwater  samples  is 
not  well  understood  at  the  present  time,  Perry,  Grundl ,  and  Gilkeson  (in 
press)  have  concluded  tentatively  that  it  is  due  to  mixing  between  ground- 
waters of  two  end  member  compositions.  This  mixing  may  be  due  to  physical 
properties  of  the  aquifer  matrix  that  affect  the  miscible  displacement  of 
groundwater  by  an  invading  groundwater  of  a  different  chemical  composition. 

The  greatest  enrichment  in  3"*S  is  in  groundwater  samples  from  well 
sites  located  in  the  strong  reducing  zone  (Kane-2  and  Kane-3,  table  1A)  and 
in  the  transition  zone  (Cook-2  and  Kane-4,  table  1A)  that  is  downgradient 
of  the  strong  reducing  zone.  Determination  of  the  fractionation  of 
<531*S/Sq  2-\  that  has  occurred  as  a  result  of  bacterial  reduction  of  sulfate 
in  the  strong  reducing  zone  requires  determination  of  the  <$31tS,S0  2-\  of 
the  initial  sulfate.  Information  on  the  isotopic  composition  of  dissolved 
sulfate  before  significant  reduction  has  occurred  is  provided  by  study  of 
the  covariation  of  <$34S/S0  2-\  and  <5180/so  2-\.  From  experiments  on  bacte- 
rial reduction  of  sulfate,  Mizutani  and  Rafter  (1969)  concluded  that,  at 
any  stage  of  reduction,  the  ratio  of  180  enrichment  to  3US  enrichment  in 
the  residual  sulfate  is  approximately  1:4.  This  relationship  is  not 
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apparent  in  the  analyses  performed  in  the  present  research.  The  <5180/Sq  2-\ 
for  groundwater  samples  from  well  sites  in  the  strong  reducing  zone  (+2.50% 
at  Kane-2  and  +4.77%0  at  Kane-3)  are  not  significantly  different  from  the 
median  value  of  +3.75%0  obtained  from  analyses  of  sul fide-free  groundwater 
samples  from  seven  well  sites  in  the  oxidized,  unconfined  region  of  the 
aquifer  system  in  north-central  Illinois.  However,  the  <5180/«.q  2-\  values 
for  the  samples  from  the  strong  reducing  zone  differ  significantly  from 
those  for  samples  from  well  sites  downgradient  of  the  strong  reducing  zone 
(+17.50%  at  Kane-4,  +14.80%,  at  Cook-2,  and  +15.91%.  at  Cook-6).  This  data 
leads  to  the  interpretation  that  the  significant  source  of  the  dissolved 
sulfate  in  waters  at  well  sites  in  the  strong  reducing  zone  is  dissolved 
sulfate  transported  into  the  reducing  zone  from  the  upgradient,  unconfined 
region  of  the  aquifer  system.  The  enrichment  in  3l*S  is  +40.23%o  for  well 
site  Kane-2  and  +29.04%o  for  well  site  Kane-3,  compared  to  the  median 
63l*S/S0  2-\   value  of  +5.07%o  for  seven  sul  fide-free  groundwater  samples  from 
well  sites  in  the  oxidized,  unconfined  region  of  the  aquifer  system.  Frac- 
tionation of  this  magnitude  may  occur  through  bacterial  reduction  of  sulfate 
(Claypool  et  al.,  1980). 

The  observation  that  no  significant  fractionation  of  <5180/<;q  2-\  has 
occurred  in  the  strong  reducing  zone  is  problematic  and  warrants  further 
study.  Fractionation  may  be  concentration-dependent,  and  the  insignificant 
isotopic  fractionation  of  oxygen  may  be  the  result  of  the  low  concentrations 
of  dissolved  sulfate  in  the  strong  reducing  zone. 

DISCUSSION 

The  isotopic  composition  of  the  water  molecule  in  groundwater  of  the 
unconfined  region  of  the  Cambrian-Ordovician  aquifer  system  is  evidence  of 
the  importance  of  glacial  processes  to  the  recharge  of  groundwater  to  the 
flow  system.  The  role  of  glacial  processes  as  a  source  of  dissolved  sulfate 
must  also  be  considered.  The  mechanics  of  recharge  related  to  glacial  pro- 
cesses are  not  clearly  understood;  the  mechanism  proposed  in  this  paper  is 
that  subglacial  recharge  of  basal  meltwater  occurred  during  the  recent 
glacial  stages  in  Illinois  (approximately  13,000  years  before  the  present) 
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when  glaciers  occupied  Lake  Michigan  with  ice  boundaries  in  the  northeast- 
ern corner  of  the  state.  As  a  result  of  the  geothermal  gradient  as  well  as 
the  heat  generated  by  friction,  meltwater  was  present  at  the  base  of  the 
continental  glaciers  (Weertman,  1972;  Clayton  and  Moran,  1974;  Freeze  and 
Cherry,  1979).  It  is  postulated  that  pressure  from  the  glacial  ice  over- 
burden caused  basal  meltwater  containing  sulfate  and  other  salts  to  re- 
charge into  the  Cambrian-Ordovician  aquifer  system  through  the  floor  of 
Lake  Michigan.  The  subglacial  recharge  reversed  regional  hydraulic  gradi- 
ents and  the  meltwater  flowed  westward  across  the  region  to  zones  of  dis- 
charge in  the  unconfined  region  of  the  aquifer  system.  Thus,  the  dissolved 
sulfate  was  derived  from  the  glacial  erosion  of  evaporites,  and  in  addition, 
from  recharge  through  evaporite  sequences  in  the  rock  strata. 

Extensive  evaporite  sequences  are  present  east  of  Illinois  in  strata 
of  the  middle  Devonian  and  the  middle  and  upper  Silurian  rocks  that  under- 
lie and  flank  Lake  Michigan  (Shaver  et  al . ,  1971).  At  present,  the  iso- 
topic  composition  of  sulfate  in  these  strata  is  unknown.  However,  analyses 
in  other  areas  of  the  world  (fig.  7)  have  established  a  &3kS(cr\   2-\  range 
of  approximately  +22%0  to  +27%0  for  rocks  of  middle  and  upper  Silurian  age. 
Analyses  of  strata  of  middle  Devonian  age  have  determined  a  <53,4S/Sq  2-\ 
range  of  approximately  +17%0  to  +19%0  (Claypool  et  al.,  1980).  Although 
sulfur  isotopic  analyses  of  evaporites  in  the  Lake  Michigan  Basin  are  needed 
to  confirm  the  postulated  mechanism,  the  available  data  strongly  indicate 
that  the  isotopic  composition  of  sulfates  in  strata  of  the  Lake  Michigan 
Basin  will  be  similar  to  the  isotopic  composition  of  sulfate  derived  from 
dissolution  of  sulfate  minerals. 

The  hypothesis  explains  several  characteristics  of  groundwater  in  the 
unconfined  region  that  is  downgradient  of  the  strong  reducing  zone  (table 
1A).  These  include  the  isotopic  composition  of  dissolved  sulfate,  the 
depletion  of  180  in  the  water  molecule,  the  high  oxidation-reduction  poten- 
tial of  the  groundwater,  and  the  low  concentrations  of  dissolved  organic 
carbon. 

The  presence  of  oxidized  groundwater  in  the  confined  region  of  the 
flow  system,  downgradient  of  the  strong  reducing  zone,  is  anomalous.  The 


-  30 


general  evolutionary  trend  of  groundwater  in  confined  flow  systems  is 
toward  a  decrease  in  the  redox  potential  along  the  direction  of  flow  to  an 
anaerobic  state;  groundwater  further  along  the  flow  path  remains  anoxic. 
However,  if  the  oxidized  groundwater  was  recharged  by  a  glacial  process 
through  the  floor  of  Lake  Michigan,  then  the  groundwater  did  not  pass 
through  the  strong  reducing  zone.  The  strong  reducing  zone  is  thought  to 
have  been  present  in  interglacial  flow  systems.  Recharge  of  groundwater 
with  a  distinct  chemical  composition  due  to  glacial  processes  may  have  had 
a  significant  effect  on  chemical  reactions  in  this  zone.  The  rapid  sub- 
glacial  recharge  of  groundwater  could  have  temporarily  restricted  the  acti- 
vity of  anaerobic  bacteria. 

During  passage  through  geologic  materials,  the  miscible  displacement 
of  fluids  with  differing  chemical  compositions  can  be  affected  by  the  physi- 
cal characteristics  of  the  matrix.  The  sandstone  matrix  of  the  aquifer 
system  can  be  partitioned  into  conducting  pores,  nonconducting-exchanging 
pores,  and  nonconducting-nonexchanging  pores.  The  characteristics  of  pores 
in  Ottawa  sand  (a  commercial  product  mined  from  one  of  the  sandstones  of 
the  Cambrian-Ordovician  aquifer  system)  were  determined  during  miscible 
displacement  experiments  conducted  by  Weeks,  Stewart,  and  Weeks  (1976). 
Only  about  2  percent  of  the  pores  were  nonexchanging  (totally  isolated). 
Forty  percent  of  the  pores  were  nonconducting  but  open  to  exchange  through 
diffusion.  The  diffusion  processes  were  thought  to  be  relatively  rapid 
because  exchange  was  essentially  complete  after  displacement  of  about  two- 
pore-volumes. 

In  a  broad  region  of  the  Cambrian-Ordovician  aquifer  system  of  north- 
eastern Illinois,  the  postulated  subglacial  recharge  would  have  resulted  in 
the  miscible  displacement  of  groundwater  with  a  high  dissolved-barium  con- 
centration by  groundwater  with  a  high  dissolved-sulfate  concentration.  The 
sulfate-bearing  groundwater  passing  through  conductive  pore  space  in  the 
aquifer  matrix  would  interact  through  diffusion  with  groundwater  containing 
dissolved  barium  in  nonconductive,  exchanging  pores.  A  widespread  condi- 
tion of  supersaturation  with  respect  to  the  solubility  product  of  barium 
sulfate  would  develop,  and  precipitation  of  finely  disseminated  barite  could 
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occur.  The  isotopic  composition  of  sulfate  in  the  barite  would  be  nearly 
identical  to  that  of  the  dissolved  sulfate,  approximately  +16%  for 
6180/S0  2-)  and  +20%  for  63ltS/S0  2-\-  Chemical  analyses  have  determined 
that  barium  occurs  widely  in  the  aquifer  strata  (Gilkeson  et  al . ,  1978). 
The  presence  of  secondary  barite  has  not  yet  been  confirmed;  however, 
sampling  groundwater  from  a  broad  region  of  the  aquifer  system  established 
that  the  activities  of  barium  and  sulfate  are  in  chemical  equilibrium  with 
the  mineral  phase  barite— evidence  that  the  mineral  is  present  in  the 
strata  (Gilkeson  et  al.,  1978). 

The  isotopic  composition  of  dissolved  sulfate  in  groundwater  samples 
from  well  sites  De  Kalb-3  and  De  Kalb-4  (table  1A  and  fig.  1)  is  evidence 
that  glacially  recharged  groundwater  was  present  to  the  limit  of  the  aqui- 
fer's confinement  by  the  Maquoketa  Shale  Group.  For  this  region  of  the 
aquifer  system,  the  isotopic  composition  of  the  low  dissolved-sulfate  con- 
centration (<10  mg/L)  is  anomalous.  The  <5180/H  Q\  of  groundwater  at  these 
well  sites  is  contemporary  with  present-day  meteoric  water.  The  primary 
source  of  dissolved  sulfate  in  these  groundwater  samples  is  interpreted  to 
be  dissolution  of  secondary  barite  precipitated  in  response  to  glacial  re- 
charge. 

Dissolution  of  secondary  barite  may  also  be  an  important  source  of 
dissolved  sulfate  in  groundwater  samples  from  well  sites  Kane-4  and  Cook-2 
that  are  located  in  the  transition  zone  downgradient  of  the  strong  reducing 
zone.  The  <53"*S/S0  2-\  composition  of  the  groundwater  at  Kane-4  and  Cook-2 
indicates  that  an  additional  source  of  dissolved  sulfate  is  3I*S-enriched 
sulfate  transported  in  groundwater  from  the  strong  reducing  zone.  Signifi- 
cantly, the  6180/™  2-\   composition  of  groundwater  samples  from  well  sites 
in  the  strong  reducing  zone  (+2.50%  at  Kane-2  and  +4.77%  at  Kane-3)  indi- 
cates that  secondary  barite  precipitated  by  glacial  recharge  is  not  present 
at  this  time.  It  is  postulated  that  bacterial  sulfate  reduction  has  con- 
sumed this  source  of  sulfate;  therefore,  a  source  of  barium  other  than 
secondary  barite  is  required  to  explain  the  high  dissolved-barium  concentra- 
tions in  groundwater  at  these  well  sites.  Depletion  of  barite  would  allow 
dissolution  of  barium  from  other  sources  in  the  strata. 
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A  factor  complicating  the  interpretation  of  analytical  data  is  that 
well  sites  Kane-2  and  Kane-3  in  the  strong  reducing  zone  are  constructed  to 
receive  groundwater  from  the  entire  Cambrian-Ordovician  aquifer  system. 
Groundwater  produced  from  these  wells  is  a  blend  from  all  the  rock  strata. 
Moreover,  these  wells  were  developed  with  explosives  to  increase  production 
of  groundwater  from  the  Ironton-Galesville  Sandstone.  Comparison  of  chemi- 
cal analyses  of  groundwater  from  well  sites  Kane-2  and  Kane-3  to  other  wells 
in  close  proximity  that  were  not  developed  with  explosives  indicates  that 
there  is  a  significant  vertical  variation  in  the  chemical  composition  of  the 
groundwater  in  the  aquifer  system;  groundwater  produced  from  the  Ironton- 
Galesville  Sandstone  has  a  lower  dissolved-sulfate  concentration  and  a  sig- 
nificantly higher  dissolved-barium  concentration  than  groundwater  produced 
from  other  strata  in  the  aquifer  system. 

At  well  sites  Kane-2  and  Kane-3,  the  composition  of  groundwater  indi- 
cates that  both  sulfate  reduction  and  methane  formation  are  occurring  in 
microenvironments  in  the  rock  strata.  It  is  significant  that  methane  is 
found  in  groundwater  at  these  well  sites.  According  to  observations  made 
in  several  studies  of  marine  sediments,  methane  does  not  occur  until  sulfate 
reduction  is  complete  (Goldhaber  and  Kaplan,  1974).  The  presence  of  methane 
is  evidence  that  sulfate  reduction  has  gone  to  completion  in  some  micro- 
environments  of  the  aquifer  system. 

After  sulfate  reduction  has  depleted  barite  (BaSCs)  in  the  rock 
strata,  another  source  of  barium  would  possibly  be  dissolution  of  witherite 
(BaC03).  Laboratory  experiments  have  determined  the  solubility  products  of 
these  mineral  phases  to  be  10"10  for  barite  and  10"8-8  for  witherite  (Sill6n, 
1964).  Witherite  has  a  greater  solubility  than  barite;  however,  in  micro- 
environments  where  sulfate  reduction  is  occurring,  witherite  may  be  less 
soluble  than  barite  because  of  the  increased  concentration  of  the  bicarbon- 
ate ion  due  to  sulfate  reduction.  In  fact,  sulfate  reduction  may  result  in 
precipitation  of  secondary  witherite.  With  the  completion  of  sulfate  re- 
duction, the  concentration  of  the  bicarbonate  ion  would  decrease,  resulting 
in  dissolution  of  secondary  witherite. 

In  an  earlier  section  of  this  paper  it  was  reported  that  confined 
groundwater  in  the  Cambrian-Ordovician  aquifer  system  in  the  west-central 
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part  of  Illinois  is  depleted  in  180.  Over  a  several -county  region  the 
6180/n  q\  ranges  from  approximately  -11%  to  -12%0 .  The  isotopic  composi- 
tion requires  that  recharge  of  the  groundwater  was  related  to  glacial  pro- 
cess. The  groundwater  in  this  region  is  highly  mineralized;  dissolved  sul- 
fate concentrations  range  from  500  mg/L  to  as  great  as  1500  mg/L.  Analyses 
of  groundwater  at  six  well  sites  in  this  region  yielded  a  median  <531*S/^q  2-\ 
value  of  +17.37%  with  a  range  from  +16.13%  to  +19.23%.  The  source  of 
the  dissolved  sulfate  is  interpreted  to  be  dissolution  of  evaporite  minerals, 
However,  similarly  to  the  situation  described  for  northeastern  Illinois, 
evaporite  minerals  are  not  known  to  occur  in  the  aquifer  strata,  and  the 
isotopic  composition  of  the  dissolved  sulfate  is  anomalous  compared  to 
values  determined  for  evaporites  of  Cambrian  and  Ordovician  age. 

It  is  proposed  that  glacial  processes  have  had  a  large  effect  on  the 
chemical  composition  of  groundwater  in  this  region  of  the  aquifer  system. 
Because  of  the  limited  information  available  at  the  present  time,  it  is 
only  possible  to  speculate  on  the  glacial  processes  involved.  During  the 
Woodfordian  Substage  of  the  Wisconsinan  glaciation,  two  glacial  events  are 
recognized  that  may  have  been  important  to  recharge  of  groundwater  into  the 
flow  system.  First,  the  Woodfordian  glacier  moved  into  Illinois  out  of 
Lake  Michigan  approximately  23,000  years  before  the  present  (Kempton  and 
Gross,  1971),  and  at  the  maximum  advance  the  ice  front  occupied  the  ancient 
drainageway  of  the  Mississippi  River  along  the  east  side  of  Stark,  Peoria, 
and  Tazewell  Counties  (Willman  and  Frye,  1970).  The  ancient  valley  (now 
occupied  in  part  by  the  Illinois  River)  may  have  been  a  significant  zone  of 
discharge  for  groundwater  in  the  Cambrian-Ordovician  aquifer  system.  Over- 
burden pressures  of  the  glacial  ice  sheet  would  have  reversed  regional  hy- 
draulic gradients  and  recharged  basal  meltwater  into  the  flow  system.  Sub- 
glacial  recharge  could  have  occurred  during  the  approximate  2,000-year 
interval  that  glacial  ice  was  present  in  this  setting. 

The  second  glacial  event  occurred  later  with  the  advance  into  north- 
central  Iowa  of  the  late  Woodfordian-age  Des  Moines  Lobe  of  the  Wisconsinan 
Laurentide  ice  sheet  about  14,000  years  before  the  present  (Kemmis,  Hal  1  - 
berg,  and  Lutenegger,  1981).  The  maximum  advance  of  the  glacial  ice  was  to 
the  vicinity  of  Des  Moines,  Iowa.  For  approximately  1,000  years,  the  ice 
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sheet  was  situated  over  the  principal  recharge  area  for  strata  of  the 
Cambrian  and  Ordovician  systems  in  Iowa  (Horick  and  Steinhilber,  1978). 
The  glacier  overlaid  bedrock  that  contains  marine  evaporite  minerals  of  mid- 
dle Devonian  and  Jurassic  age.  The  631*S  composition  of  evaporites  of  middle 
Devonian  and  Jurassic  age  (fig.  7)  are  similar  to  the  isotopic  composition 
of  dissolved  sulfate  in  groundwater  of  the  aquifer  system  in  west-central 
Illinois. 

Further  study  of  the  Cambrian-Ordovician  aquifer  system  in  western 
Illinois  and  its  relation  to  the  flow  system  in  Iowa  is  necessary  to  deter- 
mine the  role  of  glacial  processes  with  respect  to  the  composition  of  the 
groundwater.  The  multiple  glaciations  of  the  Pleistocene  are  thought  to 
have  had  large  effects.  Gilkeson  and  Cowart  (in  press)  speculate  that 
glacial  processes  are  important  to  the  distribution  of  isotopes  of  the 
uranium-238  decay  series  in  the  aquifer  system.  In  the  present  paper, 
speculation  is  restricted  to  the  effects  of  the  most  recent  glaciation  (late 
Wisconsinan).  Radiocarbon  dating  of  groundwater  may  establish  the  impor- 
tance of  Wisconsinan  glacial  events  to  the  chemical  composition  of  ground- 
water in  the  flow  system. 

SUMMARY 

High  dissolved-barium  concentrations  in  the  Cambrian-Ordovician 
aquifer  system  of  northeastern  Illinois  occur  in  groundwater  that  is  anoma- 
lously depleted  in  dissolved  sulfate.  The  dissolved-sulfate  concentration 
is  an  important  control  on  the  dissolved-barium  concentration  due  to  the 
low  solubility  of  the  mineral  phase  barite  (barium  sulfate).  The  isotopic 
and  chemical  composition  of  the  groundwater  indicate  that  the  depletion  of 
dissolved  sulfate  was  due  to  bacterial ly  mediated  reactions.  Field  measure- 
ments of  the  oxidation-reduction  potential  and  analyses  of  total  organic 
carbon  determined  that  conditions  for  anaerobic  bacterial  metabolism  exist 
in  zones  of  the  aquifer  system  that  are  depleted  in  dissolved  sulfate. 
Other  characteristics  of  groundwater  in  these  zones  that  indicate  bacterial 
sulfate  reduction  are  the  presence  of  dissolved  hydrogen  sulfide,  increased 
alkalinity,  and  enrichment  in  3l4S  in  the  residual  dissolved  sulfate. 
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Isotopic  analyses  determined  two  principal  sources  of  dissolved  sul- 
fate for  the  Cambrian-Ordovician  aquifer  system  in  northeastern  Illinois. 
The  most  significant  source  of  dissolved  sulfate  in  the  unconfined  region 
of  the  flow  system  is  the  oxidation  of  sulfide  minerals,  and  the  primary 
source  of  dissolved  sulfate  in  the  mineralized  groundwater  in  the  confined 
region  of  the  flow  system  is  from  dissolution  of  marine  evaporite  minerals. 
Evaporite  minerals  are  also  the  primary  source  for  dissolved  sulfate  in 
mineralized  groundwater  in  the  aquifer  system  in  west-central  Illinois. 
The  source  of  the  marine  evaporite  minerals  is  problematic.  Evaporites  are 
not  known  to  occur  in  the  Cambrian  and  Ordovician  strata  of  the  flow  system. 
Also  the  isotopic  analyses  determined  the  dissolved  sulfate  to  have  been 
derived  from  evaporites  that  are  not  of  Cambrian  or  Ordovician  age. 

The  depletion  of  180  in  groundwater  in  the  confined  region  of  the 
aquifer  system  is  evidence  of  the  importance  of  glacial  processes  to  re- 
charge. Recharge  related  to  glacial  processes  is  proposed  as  an  explana- 
tion of  the  anomalous  isotopic  composition  of  dissolved  sulfate.  Late 
Woodfordian  ice  sheets  occupied  Lake  Michigan  with  ice  margins  in  the  north- 
eastern corner  of  Illinois.   It  is  proposed  that  the  ice  sheet  reversed 
regional  hydraulic  gradients,  and  significant  recharge  to  the  Cambrian- 
Ordovician  aquifer  system  occurred  by  subglacial  flow  of  meltwater  through 
the  floor  of  Lake  Michigan.  A  source  of  dissolved  sulfate  is  marine  evap- 
orites of  Silurian  and  Devonian  age  that  underlie  Lake  Michigan. 

At  the  maximum  advance  of  the  Woodfordian,  glacial  ice  was  present 
over  the  ancient  Mississippi  River  Valley  in  central  Illinois.  In  this 
setting,  subglacial  recharge  of  basal  meltwater  may  have  occurred  to  the 
Cambrian-Ordovician  aquifer  system  in  west-central  Illinois.  However,  the 
most  recent  glacial  event  that  may  be  related  to  recharge  of  groundwater  to 
the  flow  system  in  west-central  Illinois  was  the  advance  of  the  Des  Moines 
Lobe  over  the  principal  region  of  recharge  into  the  Cambrian  and  Ordovician 
strata  in  northern  Iowa.  The  confined  flow  system  in  west-central  Illinois 
may  be  an  extension  of  the  flow  system  in  Iowa.  Possible  sources  of  dis- 
solved sulfate  are  evaporites  of  Jurassic  age  present  in  northern  Iowa  or 
evaporites  of  middle  Devonian  age  present  in  central  Iowa.  Further  research 
is  necessary  to  establish  the  relation  of  glacial  processes  to  the  composi- 
tion of  groundwater  in  the  flow  system. 
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